Massive triplet excitations in a magnetized anisotropic Haldane spin chain. 
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Inelastic neutron scattering experiments on the Haldane-gap quantum antiferromagnet 
Ni(C5Di4N2)2N3(PFe) are performed at mK temperatures in magnetic fields of almost twice the 
critical field H c applied perpendicular to the spin cahins. Above H c a re-opening of the spin gap 
is clearly observed. In the high-field Neel-ordered state the spectrum is dominated by three dis- 
tinct long-lived excitation branches. Several field-theoretical models are tested in a quantitative 
comparison with the experimental data. 



One-dimensional (ID) integer-spin antiferromagnets 
(AFs) are famous for having a disordered "spin liquid" 
ground state and an energy gap A ~ exp(— ttS) in the 
excitation spectrum 0. Elementary excitations are a 
triplet of massive (gapped) long-lived "magnons". An 
external magnetic field modifies the magnon energies by 
virtue of Zeeman effect 0> 0] ■ At a certain critical field 
H c the gap in one of the branches approaches zero 0, . 
The result is a condensation of magnons 0, Q and the 
emergence of a qualitatively new ground state. What 
this new ground state actually is, depends on the sym- 
metry of the problem. Theory predicts that in an axi- 
ally symmetric (AS) Heisenberg or XY-like scenarios the 
high-field phase is a gapless "Luttinger spin liquid" with 
quasi-long-range order and a diffuse continuum of exci- 
tations (no sharp magnons) The high-field phase in 
the axially asymmetric (AA) case is expected to be to- 
tally different. Here the ground state should have true 
long-range Neel order ("spin solid"). A simple boson de- 
scription Q predicts a re-opening of the gap at H > H c , 
and implies the restoration of a single-particle excitation 
spectrum. Is this state then similar to a classical easy- 
plane AF in a field, that also features long-range order 
and sharp gapped spin waves? 

Only recently did experiments, which are the key to un- 
derstanding the high-field behavior, become technically 
feasible. This was in part due to the discovery of the very 
useful model material NiCCsD^VgNs^Fe) (NDMAP) 
, where various techniques 0, El confirmed a quan- 
tum phase transition at an easily accessible critical field 
of H c « 6 T. Inelastic neutron studies were carried out 
in the AA geometry in the thermally- disordered phase: 
at H > H c , but at temperatures high enough to de- 
stroy long-range Neel order (loj . Somewhat unexpect- 
edly, it was found that as the Haldane gap closes at H c , 



the spectrum retains a considerable quasielastic (gapless) 
component at higher fields. The theoretically predicted 
reopening of the gap was thus not observed. At the time, 
this behavior was not fully understood, though several in- 
triguing explanations were put forward. One attributed 
the phenomenon to the ID diffusion of thermally excited 
classical solitons ^l), while another drew parallels with 
the incommensurate Luttinger liquid state in the AS ge- 
ometry. To better understand this issue, we carried out 
a new series of measurements at considerably lower tem- 
peratures and in higher magnetic fields, overcoming any 
finite-T effects to directly probe the ground state prop- 
erties. The spin dynamics in this regime was found to 
be qualitatively different from that previously seen at el- 
evated temperatures. The new data allow a quantitative 
comparison with several quantum field-theoretical mod- 
els, while emphasizing dramatic differences between the 
high-field phase and a classical magnet. 

Five newly-grown deuterated NDMAP single crystals 
were co-aligned by neutron diffraction to produce a sam- 
ple of total mass 1.4 g and a mosaic spread of 3°. 
NDMAP crystallizes in the orthorhombic space group 
Pnmn. The 5=1 AF spin chains, formed by Ni 2+ ions 
bridged by azido-groups, run along the crystallographic 
c axis, and (a, b) is the magnetic easy plane. In zero field 
the Haldane gap energies were previously determined to 
be A x = 0.42 meV, A y = 0.52 meV, and A z = 1.89 meV 
. In our experiments the sample was mounted with the 
a axis vertical, and the data were collected in the (0, k, I) 
reciprocal-space plane. The sample environment was a 
cryomagnet with a dilution refrigerator. The measure- 
ments were performed at T = 30 mK in magnetic fields 
of up to 11 T applied along the crystallographic a axis. 

The first series of experiments was performed at the 
SPINS 3-axis cold-neutron spectrometer installed at the 
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FIG. 1: A series of constant-^ scans measured in NDMAP 
at T — 30 mK for different values of magnetic field applied 
along the a axis (symbols). The lines are fits to a simple 
single-mode cross section function as described in the text. 



NIST Center for High Resolution Neutron Scattering 
(CHRNS). The main purpose was to measure the field de- 
pendence of scattering at the ID AF zone-center I = 0.5. 
Neutrons with a 3.7 meV fixed-final energy were used 
with a horizontally focusing analyzer and a BeO filter 
after the sample. Energy scans were performed on the 
(0, k, 0.5) reciprocal-space rod. The wave vector transfer 
perpendicular to the chains was continuously adjusted 
to maintain the sample c axis directed towards the ana- 
lyzer for optimal wave vector resolution along the chains. 
The background (featureless and typically 4 counts/min) 
was measured away from the ID AF zone-center, at 
(0,fc,0.35) and (0, A;, 0.65). 

Typical background-subtracted data sets are shown in 
Fig. ^ A similar scan previously measured in zero field 
(Fig. 3 in Ref. 0) clearly shows two peaks at roughly 
0.47 mV and 1.9 meV energy transfer, respectively. The 
data plotted in Fig. corresponds to H = 3 T, still 
well below H c 6 T. At this field the lower-energy peak 
is visibly split in two components. At H = 6 Tk H c 
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FIG. 2: Measured field dependence of the gap energies in 
NDMAP at T = 30 mK and H applied along the crystallo- 
graphic a axis (open symbols). Dashed and dash-dot lines 
are predictions of the theoretical models proposed in Refs. |3) 
and 3, respectively. The solid lines are a best fit to the data 
using an alternative model described in this work. 



the gap in the lower mode vanishes altogether, to within 
experimental error, as illustrated in Fig.^3. At the same 
time, long-range AF order sets in |8|]. The main new 
result of the present study is the observation that at T = 
30 mK at H > H c the gap in the lower mode re- opens 
and increases with field (Fig.^i,d). In the ordered state 
the spectrum at g|| = tt thus contains three distinct sharp 
excitation branches, just as in the low-field disordered 
phase. 

The new 30 mK data clearly show that the quasielastic 
scattering previously observed at H > H c and T = 2K 
is absent, and must therefore be a finite-T effect. In the 
constant-g scans collected at T = 30 mK all three peaks 
have resolution-limited widths at all fields. In fact, very 
good fits to the data (solid lines in Fig. ^) can be ob- 
tained using a simple model cross section that involves 
three excitations with a zero intrinsic energy width 
The cross section was numerically convoluted with the 
spectrometer resolution function, and the adjustable pa- 
rameters at each field were the three gap energies and 
intensity prefactors for each mode. In Fig. ^ the partial 
contributions of the three branches are shown as shaded 
areas. The field dependence of the gap energies deduced 
from these fits is shown in Fig. [3 in open symbols. 

To study the wave vector dependence of the dynamic 
structure factor additional measurements were performed 
using the Disk Chopper Spectrometer (DCS) at NIST 
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FIG. 3: Inelastic spectra measured in NDMAP using the Disk 
Chopper spectrometer for several applied fields. The range of 
the false-color scale in the lower panel is to 15 arb. u. 
Contour lines are drawn with 3 arb. u. steps in all panels. 
Arrows indicate the gap energies for the different excitation 
branches. 



CHRNS. The data were collected using a fixed incident 
neutron energy of 4.5 meV. The sample was mounted 
with the (6, c) plane horizontal, and the chain axis almost 
perpendicular to the incident beam. The background was 
measured separately, with the sample removed from the 
cryostat. The background-subtracted data collected at 
H = 0, H = 6 T and H = 10 T are visualized in the 
false-color and contour plots in Fig. |3 and correspond 
to a typical counting time of 20 hours. They are to be 
compared to similar 3-axis data measured previously at 
T = 2 K and shown in Fig. 2 of Ref. [HJ|. The new low-T 
experiment shows that the excitations at H > H c have a 
simple relativistic (hyperbolic) dispersion relation, with 
a spin wave velocity equal to that at H < H c (see solid 
lines in Fig. |3J). The "inverted" hyperbolic dispersion 
curves with "negative gaps" shown in solid lines in Fig. 3 
of Ref. ^3 are clearly inconsistent with the new data. 
This latter dispersion form was proposed as one possible 
interpretation of the anomalous g-width of quasielastic 
scattering at T — 2 K [ljj and is based on a Fermion 



representation of excitations in isotropic Haldane spin 
chains. The present data show that this model does not 
apply in highly anisotropic case of NDMAP: the spec- 
trum is truly gapped and has no intrinsic incommensu- 
rate features. The quasielastic scattering at T = 2 K, 
H > H c is thus to be attributed to a diffusion of ther- 
mally excited topological solitons and its anomalous 
q- width is due to the T-dependent mean distance between 
solitons 

For the following discussion of the observed low-T 
properties it is crucial to note that at T = 30 mK the 
spin chains are antiferromagnetically ordered at H > H c 
|8(, with a static staggered magnetization as large as 
m ~ 1/xb per site at H = 11 T. Clearly, inter-chain 
coupling is needed to stabilize order at a non-zero tem- 
perature. However, in the AA geometry, even an iso- 
lated chain orders at H > H c at T = 0, the system be- 
ing equivalent to the (l-fl)-dimensional Ising-model 0). 
Since inter-chain interactions in NDMAP are very weak 
:9j, we can assume that a purely ID problem is realized: 
long-range AF correlations are intrinsic to the ID chains, 
and the sole role of residual 3D interactions is to maintain 
their stability at a finite temperature. 

The conventional approach to describing spin excita- 
tions in ordered systems is the quasiclassical spin wave 
theory (SWT). In this model the magnons are precessions 
of staggered magnetization m around its equilibrium di- 
rection. As a consequence, in SWT there are only two 
sharp excitation branches, polarized perpendicular to the 
m. In our case three sharp magnons are seen in the or- 
dered state (H > H c ). At least one of the three branches 
must have the character of a "longitudinal" magnon that 
is not a precession mode, but is polarized along the or- 
dered moment. Thus, at H > H c quantum-mechanical 
effects remain crucial, and the SWT is inapplicable. In- 
stead, the three observed excitation branches can be visu- 
alized as soliton-antisoliton breathers: the three massive 
bound states formed by the two types of topological de- 
fects allowed in an anisotropic semisiclassical ID magnet 
|ll| . In more detail our experiments can be understood 
in the frameworks of several field-theoretical models. 

The approach due to Affleck j^] is based on coarse- 
graining the (l+l)-dimensional 0(3) non-linear sigma 
model (NLSM), to which, in turn, the 5 = 1 Heisen- 
berg chain can be mapped Q . The resulting Lagrangian 
is that of an unconstrained real vector field ip(r) with the 
<p -type interaction. Anisotropy is added by postulating 
separate masses A Q for the different components of this 
vector field. For H > H c the ground state has a non-zero 
staggered magnetization L = (tp) and uniform magneti- 
zation M cx (H X <p). The iy9 4 -model captures the ba- 
sic physics involved, but suffers from several drawbacks. 
In particular, the predicted value of the critical field for 
H||e a is gusHc = A „. T his is inconsistent with estab- 
lished experimental @,13E3>E1E| an d numerical 0] 
results, as well as with simple arguments based on the 
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perturbation theory [lTj. Another potential weak- 
ness is that at the mean field (MF) level M(H) and the 

lowest gap A(H) come out to be oc \j H — Hc a \ while 
a roughly linear behavior is seen experimentally and nu- 
merically. For NDMAP, the derivations of Ref.y must be 
generalized to allow for an arbitrary field direction, since 
the anisotropy axis of the Ni 2+ ions forms an angle of 
about 16° with the crystallographic c axis. The resulting 
predictions are shown in dashed lines in Fig. [21 

Another theory, due to Tsvelik Q, stems from the in- 
tegrable SU (3) model of a 5 = 1 chain. It involves three 
Majorana fields with masses A a . The predicted critical 
fields g^BH^ = \/ A g A^ coincide with the perturba- 
tive formulas of [HI 113 • Again, incorporating an arbi- 
trary field direction, and using the known parameters for 
NDMAP, we can directly compare the predicted gap en- 
ergies (dash-dot lines in Fig.[5J) to those measured in this 
work. Unlike the ip 4 model, Tsvelik's approach predicts 
the linear behavior of M(H) and A(H) near H c . The 
theory, however, fails to reproduce the two upper gaps at 
H>H C . 

In the context of the present experiments we would 
like to introduce a different approach that is similar to 
the model proposed for dimerized S — 1/2 chains in 
Ref. It is a Ginzburg-Landau-type theory written in 
terms of a complex triplet field 3>(r) = A(r) + iB(r). 
The uniform and staggered magnetization are written 
as M cx (A x B) and L cx A(l - A 2 - B 2 ) 1 ' 2 . In- 
teraction terms (3>* ■ <I>) 2 and ($* x 3>) 2 are natu- 
rally present in the model. Unlike the models discussed 
above, anisotropy enters the Lagrangian through two 
sets of masses, separately in the A and B "channels": 
J2 a ( ma ^a + ">TiaB 2 ). By integrating out the B-ficld 
one obtains an effective (p 4 -type theory similar to that 
of Affleck. However, the Zeeman term H ■ (ip x dtf) 
of Ref. becomes replaced by the anisotropy-dependent 
expression ^ a ^ 7 (rn^) _1 e Qi a 7 i? Q A^(i9 t ^4 7 ). Remarkably, 
this model includes Affleck's theory as the special case 
rrCy = const, and at the same time it reproduces the re- 
sults of the Tsvelik's theory for A(H) and M(H) below 
H c in another special case = m 7 . There are no par- 
ticular reasons why either of these special cases should 
correspond to the actual Heisenberg 5=1 chain with 
single-ion anisotropy realized in NDMAP. The masses 
and m a at the present stage may be treated as ad- 
justable parameters to reproduce the measured gap en- 
ergies A a — y/rn^rria at H = and the measured field 
dependencies. Very good fits to our experimental data 
are obtained with m x /m x = 0.87, rn y /m y = 0.83, and 
rn z /m z — 0.35 (solid lines in Fig. EJ. A detailed descrip- 
tion of the application of this model to NDMAP will be 
reported elsewhere. 

In summary, the present low-T study is a direct obser- 
vation of several fundamental quantum-mechanical fea- 
tures predicted for anisotropic Haldane spin chains, and 



helps clarify the nature of the previously investigated 
finite-T spin dynamics. The next experimental challenge 
will be to investigate the AS geometry, searching for man- 
ifestations of the Luttinger spin liquid regime. 
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